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Abstract 

The efficacy of the modelling of the solidification of castings is typically tested against 
observed cooling curves and the final grain structures and sizes. Without thermo-solutaJ con- 
vection, cquiaxed grain formation is promoted by introduction of heterogeneous substrates 
into the melt, as grain refiners. With efficient thermo solutal convection, dendrite fragments 
from the mushy zone can act as an intrinsic source of cquiaxed grains and resort to grain 
refining additions is unnecessary. The mechanisms of dendrite fragmentation and transport of 
these fragments are briefly considered. 
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INTRODUCTION 

s=i=^sr-£S 

PHiSSi 
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becomes possible, we shall need to get the physics right and it seems apptopnate to examine 
some aspects of this many sided problem, with that in mind. 

MODELLING PROCEDURES 

Se elements whfch, although cooling or locally heating, as during rccatescence, ^ 
regarded as being spatially isothermal, and, by implication, with 

(Chennai di tTus.v^ty^K, wl (iniw ckrMnl^r^OO^un^whicto may be 

inconsistent tacets oi ui v raJC of heat rem0 val is constant, or can be cal- 

follows that originally demonstra ed by Uiane ; nucleation rate. 

dence becau lhe source ( s ) 6 of nuclei and their activity(ies) vs undercoolings are not 

weEown and this is of concern here. Either, grains nucleate on heterogeneous substrates 
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.0 predict the microstructure of a fab- 
on to predict optimal mechanical or 
e best processing path. To be able to 
\ intellectually satisfying and it obvi- 
e. In the metallurgical industry, cast- 
r path for producing final shapes or 

0 the challenge is, to be able to pre- 
this stage. 

wth, sometimes taking place concur- 
d by the rate of heat removal, prima- 
e material, a metal, the most direct 
ification reaction is that of latent heat 
rse, by examination of the final grain 

1 the overall process must be able to 
nd the accuracy of this curve fitting, 
res is the measure of success. We are 
ti actual cooling curves to adjust the 
d in many cases (eg. 1-3) the results 

first principles, if indeed that ever 
and it seems appropriate to examine 
lind. 

URES 

jt the distances over which this takes 
on points. In the following, attention 
h determine the final grain size in the 
:ture would be one of fine, randomly 
i, with attendant problems of macro- 
ed. The heat flow analysis, although 
' linking finite elements of “suitable” 
jrmal analysis results also relate to 
eating, as during recalescence, are 
ation, of dimensions less than Jxt 
r a liquid metal (k * lO^mV 1 ), with 
element of K 300|im, which may be 
)n a statistical basis, these physically 
probably be overlooked and in the 
at removal is constant, or can be cal- 
:nt environment. The procedure then 
the nucleation and radial growth of 
ply assumed a given nucleation rate, 
*i, implicitly assuming that they were 
iron or graphite fragments) which 
these could be predicted with confi- 
^ very well documented in cast irons. 
ctivity(ies) vs undercoolings are not 
lucleate on heterogeneous substrates 



(a) 


Fig. 1(a) Initial part of cooling curve, 
computed for time dependent 
nucleation of aluminium on a 
substrate dispersion (nomi- 
nally ofAl 3 T 1# with 0 = 7°) of 
particle density, N p /v = 
10 3 mm'^ and radii, r = l|im. 



Fig. 1(b) Enlarged region for recales- 
cence, showing superim- 
posed nucleation rate, N, and 
sum of nucleation events to 
give a final grain size, N G /v. 


2 4 20 32 36 

(b) TIME, s 



tide, showing excess of substrates for different assumed 
contact angles, 0. Refer to Fig. 1 for 0 = 7°. 


567 





which are deliberately introduced to the melt, as grain refining additions, as in alummium 
base castings (eg. 5 , 6 ), or, in the absence of such identifiable sources, the most probable ori- 
j gins are taken to be particles of the base material, intrinsically present in the system as crys- 

tal fragments, especially, in larger castings, dendritic fragments from an earlier developing 

mushy region. 

i 


In the first case, the nucleation occurs only while the temperature is falling and a little after 
the recalescence has taken place, the grain size is determined at this very early stage and 
there may remain a considerable excess of particles which did not have an opportunity to act 
as substrates. Figs. 1 and 2, summarize this situation (7), although the action of the substrates 
remains obscure and difficult lo quantity (eg. 6 - 10 ). 

■ 

In the second case, with a dispersion of surviving primary fragments, there is no nucleation 
in the classical sense and the recalescence is entirely described by growth. It follows, that all 
the surviving crystal fragments grow into grains and there is no residue. If, as it is generally 
assumed and is well documented (eg. 11-16), these fragments originate from within a pro- 
ceeding columnar mushy zone, it is necessary to address the problem of (i) how do the frag- 
ments develop within the mushy zone and (ii), how are they transported from the mushy 
region into the bulk liquid? These two aspects of the intrinsic model will be considered here. 


DENDRITE FRAGMENTATION 

In a casting, following the evolution of a columnar front from an initial chill zone, the den- 
dritic growth proceeds as a well oriented array (< 100 > in cubic metals), with the growth rate 
decreasing (approximately, parabolically) and the temperature gradient falling, as the thick- 
ness of the casting develops. The growth conditions are therefore always transient. 

■ 

! The geometry of the growth, Fig. 3, is described by a primary interdendritic spacing, X h 

f measured across a square grid section and a primary dendrite tip radius, r lt both dimensions 

determined by the growth rate and temperature gradient, leading to an undercooling, AT, at 
\ the growth front, for a given concentration in a given system (eg. 17, 18). The undercooling is 

|j almost entirely produced by solute build up, AC, and the curvature undercooling is negligi- 

1 f bte. A stable, steady state solution to these dimensions is preferred. 

Behind the growth front, secondary dendrite side arms devclope with a nascent spacing, X°, 

! • which, for a given material, is in an approximately constant ratio to rj (19). For example, in 

; \ NH 4 CI-H 2 O, the ratio is * 4X. Ripening of the side arms begins, almost at once, to give a 

coarser secondary spacing, according, approximately, to a relationship of the form (k l 2 ) 

jh . (X°) 3 = Cjt (eg. 20 ), where Ci is a constant and t = time. 

1 1 

T Under steady state conditions, without interdendritic fluid flow. Fig. 4, the ripening assembly 

h shows that surviving side arms typically develope a narrow neck, where they join the pri- 

jj mary stems, but remarkably few, if any, actually become detached. This observation is, at 

I first sight, perplexing, because the evidence is that very large numbers of dendrite fragments 
can be seen to escape from the mushy region of a casting (eg. 21). It is therefore necessary to 

I I address this contradiction. 

[ 

■ 

The immediate reaction, knowing that there can be (see next section) vigorous reentrainment 
and circulation of bulk liquid into the mushy region, is to suppose that this is the cause of 
i enhanced side arm detachment, either by mechanical deformation at the narrow necked posi- 

j tions, or by accelerated melting, because the bulk liquid is warmer (positive temperature gra- 
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Fig. 3 Schematic dendritic growth front to show relevant dimen- 
sions which are referred to in text. 



Fig. 4 Dendrite growth in NH 4 CI-H 2 O 
at 20 pms _1 under steady state 
conditions - without side arm 
detachment 



Fig. 5 Similar region, growth at same rate, but 
followed by deceleration to 8 jims 1 
- showing extensive detachment of side 
arms. 
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I 

dient). 

i 

Of these two possible mechanisms, the first can probably be discounted, certainly for metals, 
because at the melting point, the material is so ductile that arms could only bend, but not 
fracture. We have no evidence of bent side arms. 


; The second possibility, that melting off is caused by entrainment of the warmer liquid, is 

| more attractive, until it is remembered that the bulk liquid is also less concentrated than the 

infprHpndriric liquid, i.e.. it is solvent rich, so that as it cools, solid will be deposited (prefer- 
entially into the flow direction) and melting will not be accelerated. This then is not an obvi- 
ous cause of the side arms becoming detached. 

1 

1 

i , 

The solution to this problem must lie in the fact that all real growth conditions are transient - 
as noted, with growth rate, V and temperature gradient, G, falling continuously. It is neces- 
sary to examine how the growth front, and the region behind it, respond to these changes. 
Observation shows, Fig. 5, that with a step reduction in growth rate, without convection, a 
i majority of side arms become detached - almost at once. Jackson et al. showed this effect 

( nearly twenty years ago (12), although they did not pursue it further at that time. 

* 

< 

1 

s 1 
;l 

, i 

■ f ■ 
1 ' 

, i 

[ 

The observation then is, that at the growth front, the radius of curvature, r^ responds imme- 
diately to change in growth rate, increasing as velocity decreases, according to Vr^ = con- 
stant In contrast however, the primary spacing, Xj, does not change rapidly, because it is 
constrained by the orthogonal restrictions which are created by the small anisotropy of the 
interfacial energy. does not increase immediately, with rj, nor, with acceleration, does it 

decrease rapidly. To increase Xj, it is first necessary that some dendrites fall back behind the 
front and that others adjust by subsequent branching. To decrease Xj, new primary dendrites 
can only be created from the array by development of tertiary dendrite arms. There is conse- 
quently an inherent hysteresis to adjustments of the primary spacing. Fig. 6, from Jackson et 
al. shows this very clearly. A’ - velocity increase, B' - decrease to the original. 

‘ : 

■ ; i 

i t 

r i 
\ 

How, exactly, the reduced velocity causes side arm detachment is complex. The remelting 
can be caused by (a) a rise in temperature and (b), a rise in concentration of the interdendritic 
liquid. Analysis of array growth (18) under steady stale conditions, indicates that the pre- 
ferred stable conditions are when there is a limited interaction of solute rejection from adja- 
cent dendrites (k \ , <«), such that the undercooling is not at a minimum. A small reduction in 
velocity at constant spacing, is attended by a small reduction in tip undercooling (i.e. a tem- 
perature rise at the front), although the instant increase in tip radius, r t , must simultaneously 
increase solute concentration between the dendrites, which are then “too close” for preferred 
stable array growth. 

' i i 
! 

' i 
: t 

Therefore, although the analysis is complicated, this has to contain the explanation for there 
being so many detached dendrite fragments available in a typical columnar mushy region. 
“Typical”, in the sense that the growth conditions are always transient, with V and G 
decreasing. 

-1 

ii 

j 

i 

As noted elsewhere (16), the number of side arms in a dendritic array is enormous, eg. 10 6 - 
10 7 cm" 3 , so that even if only a minor fraction were to become detached, the number of frag- 
ments, potentially available, is more than adequate to promote intrinsic nucleation - if they 
can escape from the mushy zone. 

j ■: i 

' • . !l 

i 

t ; 

The effects of liquid reentrainment by convection are probably of minor or negligible impor- 
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Fig. 6 Showing dendritic array growth in CBr^ from Jackson 
et al (12), with acceleration at level A' and deceleration 
at level B', illustrating the rapid changes in tip radius, t\ % 
slower adjustment of primary spacing, Xj, and detach- 
ment following deceleration. 


C 0 Jp/dl*-(>CuiT/dl 





(a) (b) (C) (d) 


Fig. 7 The situation leading to thermo-solutal convection. 
Schematic plots for alloy of concentration, C 0 , on phase 
diagram (a), with dendrite tip concentration, Q, grow- 
ing at an array front, (b), where the concentration vs 
height along axes A and B are as in (c) and the vertical 
liquid density profile is as in (d). 
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tance in this respect, but are ail important for subsequent transport of the fragments, once 
they have become detached. 

FRAGMENT TRANSPORT - HEAT PUMPS 

The concern is with how detached dendrite arms can be transported from the mushy region 
of a casting into the bulk open liquid, where survivors may grow into new grains. 

This is possible if interdendritic liquid, richer in solute, is less dense than the bulk liquid, 
despite being cooler, and occurs by thermo- solutal convection. The situation is most clearly 
demonstrated by solidification, vertically upwards, in a base chilled configuration, and 
becomes probable when the solutal density coefficient, p, is negative, and pAC exceeds aAT, 
where AC is the composition difference between an element of interdendritic liquid and the 
bulk, a is the coefficient of thermal expansion and AT the temperature difference corre- 
sponding to AC. From the phase diagram, AT/AC = m L , the slope of the liquidus boundary. 

There is then density inversion, given by Ap = p(pAC + aAT), between the interdendritic liq- 
uid at any position, with respect to the bulk liquid, and also between the solute boundary 
layer at the dendritic front and the latter. Fig. 7. Two convective regimes dcvclope. First, at 
very short times, c.g. < 10s, short range (< 10mm) oscillatory patterns form, close to the 
growth front. With growth vertically upwards, these create a mixing zone, immediately 
above the growth front, but the interdendritic liquid is not disturbed to any depth greater than 
that of the primary interdendritic spacing, (typically * 0.2mm). At longer times, 
eg. > I0 3 s, longer range channel plumes develope, ol diameters ® I mm (3-5 primary den- 
drite spacings) and rise far (> 100mm) into the supernatant open liquid, Fig. 8. These arc fed 
by lateral reentrainment of the warmer bulk liquid, which is drawn down through the den- 
dritic mesh, cooled and recirculated into a wider interdendritic channel. Such channels are 
subsequently choked by smaller crystals and give rise to "freckel” defects in directionally 
solidified ingots (22). The channels arc also richer in solute and are recognized in steel billet 
castings as “A” segregates. 

It follows, in all of this, that the prime requirement for such developments depends upon 
there being a sufficiently large, negative solutal density coefficient, p, so that this form of 
thcrmosolutal convection is found in steels with carbon, silicon etc., or nickel base alloys 
containing aluminium or titanium, but is not to be expected in aluminium or magnesium 
alloys which are solidifying upwards, in which a majority of solutes (eg copper) are the more 
dense. This is why intrinsic nucleation is not found in such systems, as in D.C. castings, and 
why it is then necessary to use deliberate additions of heterogeneous substrates to promote 
equiaxed grains in foundry practice with these light metals. 

Returning to the channel convection, it will be evident that the region in the mushy zone, 
below a convection plume is responsible for the near steady-state recirculation and behaves 
like a pump. The density inversion follows from the phase transformation which is driven by 
heat abstraction, so that the system behaves like a local heat pump. The relatively long range 
plumes also recirculate heat and solute into the open liquid, assisting in the heat transfer and 
causing macro segregation. Typically, the rapid local convection translates about 1-2% of the 
liquid volume upwards while the balance sinks slowly. The macrosegregation above the den- 
dritic front is described by an effective partition coefficient, k^ >0.9 (23). 

To analyze such a channel plume, as in Fig. 8, a simplified model can be used, supposing that 
a thin walled tube is inserted to some depth into the mushy region, h [> and extends upwards 

572 


ransport of the fragments, once 


PUMPS 

nsported from the mushy region 
;row into new grains. 

less dense than the bulk liquid, 
on. The situation is most clearly 
base chilled configuration, and 
negative, and pAC exceeds cxAT, 
t of interdendritic liquid and the 
le temperature difference cone- 
slope of the liquidus boundary. 

/), between the interdendritic liq- 
so between the solute boundary 
ctive regimes devclope. First, at 
ttory patterns form, close to the 
ite a mixing zone, immediately 
sturbcd to any depth greater than 
y m 0.2mm). At longer times, 
leters - 1 mm (3-5 primary den- 
jpen liquid. Fig. 8. These are fed 
is drawn down through the den- 
ritic channel. Such channels are 
frecker defects in directionally 
and are recognized in steel billet 


jc h developments depends upon 
efficient, p, so that this form of 
ilicon etc., or nickel base alloys 
ted in aluminium or magnesium 
f solutes (eg copper) are the more 
systems, as in D.C. castings, and 
rogeneous substrates to promote 


it the region in the mushy zone, 
y-state recirculation and behaves 
ransformation which is driven by 
l pump. The relatively long range 
. assisting in the heat transfer and 
:tion translates about 1-2% of the 
macrosegregation above the den- 
■ kcff >0.9(23). 

nodel can be used, supposing that 
region, hj> and extends upwards 


into the open liquid, to some height h L . The resulting recirculation becomes a fountain - a 
salt fountain (24) - driven by the pressure difference, P = Ap gh D , at its base. 

Flow is essentially streamlined, in the channel-plume pipe and also through the mesh, down- 
wards - estimates of Reynolds numbers confirm this (25). The volume of liquid flowing 
upwards, V T = JiP.r 4 /8h L Ti, must equal that flowing downwards, V i, to supply liquid 
around the edge of the base of the pipe, of circumference 2rcr. Description of the flow chan- 
nel paths, through the dendritic mesh is a more complicated problem, depending on its per- 
meability and hence, upon the primary and secondary dendrite arm spacings, X\ and X 2 , and 
the fraction of solid, f s . The pressure difference, P. at the base of the pipe is the same for 
upward and downward flow. Therefore, for a given system and set of growth conditions, it 
must follow that 


C, r 4 T=C 2 ri 

where Ci and C 2 are constants and r = (C 2 /Cj) 1/3 . 

The viscosity is essentially constant, so that for any r, the volume flow rate through the foun- 
tain is dictated by the permeability of the mushy region down to the depth, h D . This, in turn, 
is dependent on the primary and secondary dendrite spacing, Xj, which is itself related to D L / 
V d , the solute diffusion coefficient and dendritic growth rate. Since liquid diffusion rates are 
generally well within an order of magnitude (around 10' 9 m 2 s‘ l ) it becomes understandable 
why channel plume dimensions (although not the flow rales in them) arc so similar from sys- 
tem to system, metallic, aqueous, organic. 

It is possible to make artificial salt fountains of this type, where the radius, r, and lengths, (h^ 
+ h D ), are predetermined, but in the unconstrained natural setting, the model is relaxed and 
the situation becomes more involved. 


The basic balance remains, i.e. the upward and downward flows must remain equal, so that r 
is still of the form, ^/Cj) 173 , but solutions to C! and C 2 are more difficult: We suppose that 
the effective pressure term remains the same for upward and downward flows, but the walls 
of the channel within the mushy region arc now permeable and receive recirculating liquid 
over a depth range, so that hp is no longer precisely defined, - there is a continuum of possi- 
ble entrainment paths, through the dendritic array. Al the same lime, there is no stationary 
rigid wall to the plume, pipe, because the surrounding liquid is drawn up with it. The essen- 
tial requirement for this type of flow is that although thermal equilibrium is approached rap- 
idly (between the cooler liquid of the mushy zone and that of the supernatant bulk), solute 
diffusion is so much slower, that the plume composition remains essentially constant for long 
limes (and therefore distances). The composition profile across a plume is therefore of a 
square wave form, while the velocity profile is smeared out over several plume diameters. 
Such a profile can be described by spline fitting a parabolic velocity profile for the plume 
core, with an exponential profile for the surroundings - Fig. (9). Analysis of such plume flow 
shows that in transparent systems (aqueous, organic) the observed flow rates can be related 
to the buoyancy forces satisfactorily, and, by extrapolation, can be used to predict plume 
flow rates in opaque metals. Although plume dimensions are very similar cf. above, flow 
rates within them arc in a ratio, metallic: aqueous: organic * 100:10:1 mms , reflecting, pri- 
marily, the respective thermal conductivities and kinematic viscosities (22) as in t he P6cklet 
numbers. The above mentioned “constant”, C h is then to some extent numerically soluble 
and supported by experimental measurement, but the reentrainment is the more difficult to 
describe. Part of this latter problem arises because the fraction of solid in the upper pan of 
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Fig. 8. Showing convection plume above dendritic growth front of 
NH4CI-H2O, from (22), solute contrast of refractive index, plume 
width < l mm. 


V 



Fig. 9 Schematic velocity profile, V, and step concentration difference, 
AC, across a plume as in Fig. 8. 
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the mushy region is quite small (eg < 10%), so that fluid flow through the mesh does not lake 
place along identifiable, discrete channels, but rather around a very permeable array of pri- 
mary and secondary cylindrical obstacles (somewhat similar to water flowing through a reed 
bed). The shortest paths for reentrainment are obviously near the channel mouth, which is 
where the cone of primary dendrites surround the cooler plume liquid. At some lower 
depth, < 5 - 10mm, the interdcndritic liquid is essentially quiescent, and the width of the 
pump activity is not more then about twice this depth - hence the typical plume spacing. 

The result of these relaxations (on the rigid pipe model) is to increase the flow rate tn a chan- 
nel plume - by making the base of the pipe more open and by removing the constraint ol 
rigid walls. Observation and calculations for the NH4CI-H2O system, indicate an increase in 
flow rate by a factor of about x 5 (see appendix B, ref. 22). Description of the interdendntic 
flow pattern and rates remains empirical, at present, and this complicates analysis of side arm 
transport from the mushy zone sec, e.g. Poirier et al (26, 27) and Beckcrmann et al (28-30). 

SUMMARY - CONCLUSIONS 

• Without thermosolutal convention to transport dendrite iragmenis from the mushy zone 
of a casting into open liquid, nuclealion of equiaxed grains is generally promoted by 
addition of grain refining substrates. The grain size can then be modelled by lime depen- 
dent, heterogeneous nuclealion, although description ol the substrate efficiency remains 
obscure and requires arbitrary assumptions. 

• With density inversion in the columnar mushy zone, thermo-solutal convection provides 
a mechanism to disperse dendrite fragments into the melt and equiaxed grain formation/ 
nuclealion can be intrinsic, subject to fragment survival in the melt 

• The detachment of dendrite side arms in the mushy region is accelerated by temperature 
rise and by solute enrichment of interdcndritic liquid. This situation is promoted by the 
typically decreasing growth rates and temperature gradients at the columnar growth front 
of a casting, and occurs because the primary dendrite spacing does not automatically 
adjust to the continuously changing conditions, while the dendrite tip radius does so, 
immediately and continuously. 

• In a vertical configuration, finger and plume convention patterns develop if solute is suf- 
ficiently less dense than solvent, and thus offers a mechanism for efficient dispersion of 
crystal fragments into the melt. This mechanism is not generally available in light met- 
als, especially in vertical configurations. 

• Analysis of plume convection allows flow rates and dimensions to be predicted with 
some confidence, describing the mechanism in terms of a heat pump. The most intracta- 
ble part of such modelling remains in analysis of intcrdendritic flow patterns. 
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